Heat capacity of Schottky type in low-dimensional spin system 
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The heat capacity of low-dimensional spin systems is studied using theoretical and numerical 
techniques. Keeping only two energy states, the system is mapped onto the two -level-system (TLS) 
model. Using the low temperature Lanczos method, it is confirmed that the behavior of Tm and the 
energy gap as functions of the control parameter is the same in the two models studied; a conclusion 
that can probably be extrapolated to the general case of any system that possesses an energy gap. 
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I. INTRODUCTION 

Gaps in the energy spectrum play a crucial role in con- 
densed matter physics. Fundamental properties in su- 
perconductivity or in the fractional quantum Hall effect 
originate from the existence of a gap between the ground 
state and the excited states. In particular, low dimen- 
sional quantum spin systems are extremely interesting to 
study the behavior of the gap. Many exact and numeri- 
cal results on the one dimensional quantum spin systems 
with nearest neighbor couplings have been accumulated 
during last decades. 

The ID spin- 1/2 system has been solved by BetheA 
in 1931 with his famous ansatz. The ansatz allows the 
computation of the energy eigenvalues. The isotropic ID 
spin- 1/2 system with nearest neighbor couplings is gap- 
less. The anisotropic spin- 1/2 chain is denoted by the 
XXZ model 2 -! 3 -. The Hamiltonian of the XXZ model on a 
periodic chain of N sites is 



N 

H = J 2^ SfSf+i 

i=l 



(1) 



where J > is the exchange coupling in the xy easy 
plane, A is the anisotropy in the z direction. The Ising 
regime is governed by A > 1 and there is a gap in the 
excitation spectrum, while for A < — 1, the ground state 
is in the ferromagnetic phase and there is a gap over the 
ferromagnetic state. In the region — 1 < A < 1, the 
ground state of the system, is in the gapless spin-fluid 
phase. 

The spin-1 system is not solvable with the Bcthc ansatz 
or similar techniques. The anisotropic spin-1 chain is 
only gapless 4 - in the region — 1 < A < 0. Haldane^ for- 
mulated in 1983 his famous conjecture that quantum spin 
chains (isotropic) with integer spin S = 1,2,.. have a gap, 
where as chains with half-integer spin S = 1/2, 3/2, ... are 
gapless. 

During the last two decades ladder-systems^ as quan- 
tum spin systems between one and two dimensions have 
been studied as well with numerical methods. Con- 
cerning the ground state properties has been found 
that, ladder-systems with an even number of legs (I = 



2,4,6,...) have a gap; those with an odd number 
(I = 1,3,5,...) do not. In particular, since the an- 
tiferromagnetic two-leg ladder systems have a gap in 
the spin excitation spectrum, they reveal extremely 
rich quantum behavior in the presence of a magnetic 
field^. Such quantum phase transitions in spin systems 
with gapped excitation spectrum were indeed studied 
experimentall y 8 ! 9 ! 10 ! 11 ! 12 ! 13 . 

On the other hand, investigating the behavior of the 
energy gap of spin systems in vicinity of quantum critical 
points has attracted much interest recentl y 3 ' 14 ! 15 ' 16 . In 
general, the critical point of an thermodynamic system 
in the Hamiltonian formulation is defined as the value at 
which the energy gap vanishes as a power low, which is 
known as the scaling behavior. The opening of the energy 
gap in the vicinity of the quantum critical point is found 
to scale with a critical exponent. The value of the critical 
gap exponent is very important to find the universality 
class of a continuous quantum phase transition. 

The discovery of gapless or gaped excitations have 
led to the investigation of the thermodynamic proper- 
ties. One of the most important thermodynamic func- 
tions is the heat capacity. Usually, there is a lambda-type 
anomaly in figure of the heat capacity versus tempera- 
ture. It was interpreted as indicating a phase transition 
to a magnetically ordered phase. An important charac- 
teristic of the low-dimensional magnets is the absence of 
the long range order in models with a continuous sym- 
metry at any finite temperature 1 ^. There is also a broad 
maximum in the plot of the heat capacity vs tempera- 
ture, characteristic of low dimensional systems. This is 
known as the Schottky peak. 

In this paper, theoretical and numerical results are re- 
ported for the low-temperature behavior of the heat ca- 
pacity in low dimensional spin systems. Theoretically, 
by keeping only two lowest states, the system is mapped 
to the well known two-level-system (TLS) model. In this 
case, the heat capacity is found exactly as a function 
of the energy gap and the temperature. It is shown that 
the position of the Schottky heat capacity peak, Tm , and 
the energy gap behaves in the same way as a function of 
the control parameter. In Section II, the mapping to the 
TLS model is explained and the theoretical results are 
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presented. In Section III, the results of the low temper- 
ature Lanczos calculations are presented. Numerically, 
Tm as a function of the control parameter is computed 
for "the alternating spin- 1/2 chains in a magnetic field 
h" and "the ID Heisenberg Hamiltonian with a staggered 
magnetic field h s " ■ Finally, the summary and conclusions 
are presented in Section IV 



II. LOW TEMPERATURE LIMIT HEAT 
CAPACITY: TWO LEVEL SYSTEM APPROACH 

In this section we discuss a theoretical approach to find 
the effect of the energy gap on the heat capacity of the 
quantum spin systems. Spectrum energy of a quantum 
system may be gapful or gapless. Since, we are going to 
study the sign of the gap on the heat capacity, gapful 
systems are considered. At very low temperature we can 
consider only lowest energy levels. If we keep only two 
lowest energy states, the system maps to the two level 
system (TLS) modeli&. We have assumed that the energy 
gap of this TLS is g where E\ = Eq + g, E and E\ are 
the ground and first exited states respectively. 

Our purpose is to determine the behavior of the heat 
capacity for the quantum spin systems. The heat capac- 
ity is expressed by the following relation 



C v — 



1 



K B T 2 



d 2 InZ 



dp 2 



where Z is the partition function denoted by 
Z = Trje-^}, 



(2) 



(3) 



Here (3 = K T , Kb is the Boltzmann constant and T is 
temperature. We have assumed that Kb = 1. At very 
low temperature limit (TLS limit), we can write 



-0E a 



_ e -0E o 



(1 
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And finally from the above equation one can shows 

,.2 



Cy 



cosh 2 (x) ' 



(4) 



(5) 



where we defined x 



2T • 



This result predict a Schot- 



tky like peak of the heat capacity behavior versus the x 
variation. The position of the Schottky peak takes place 
at xm — 1-2, where we have: tanh(a;jvf) = 1/- t m- This 
result shows explicitly an upward increase of the heat 
capacity versus the magnetic field for x < xm and mono- 
tonic decrease for x > xm- We have plotted the thermal 
behavior of Eq. [5] in Fig. [TJ Here, xm corresponds to 
the position of the Schottky peak, Tm, in a constant gap 
value. It is clear that by increasing the gap value, \g\, 
the peak approach to higher temperature and inversely. 
Therefore we can conclude that width of the energy gap 
may affect the position of the Schottky heat capacity 
peak. At very low temperature regime, all gapped ther- 
modynamic systems (N — > oo, N = number of the spins) 
can be mapped to the above TLS model. 




Temperature 

FIG. 1: Temperature dependence of the heat capacity of TLS 
(Eq. [SJ) for different energy gap values (g). Tm shows the 
position of the Schottky heat capacity peak for the system 
with g = 2. 



Up to now we did not consider any degeneracy. In the 
general case both first two energy levels (TLS) have a 
degeneracy. It is well known for two level system (TLS) 
that the account of degeneracy leads to change the effec- 
tive gap and consequently to change the Tm for Schottky 
anomaly of heat capacity 19 . In peresent of degeneracy by 
a little of manipulation, one can shows 



Cy — 



cosh 2 (a; + x') ' 



(6) 



which is very similar to Eq. [5J Here we denoted x' — — 



IT ' 



and g s = T In ^ . Where d\ and di have been con- 
sidered for the order of the degeneracy of the ground 
state and first exited state of the system respectively. 
As we mentioned it before this result predict a Schot- 
tky like peak in the heat capacity behavior. The Schot- 
tky heat capacity peak takes place at x m which satisfied 
tanh(a;M + x' M ) = 1/xm, where x' M — 27^7- If we as ~ 
sumed that d± < di then g a < and therefore the gener- 
alized gap (g + g s ) will be smaller than the gap. Thus the 
Schottky peak moves to the higher values of the gap in 
respect to the case g s — (without degeneracy) . In same 
manner for the case: di < d\ we have g s > 0, therefore 
the generalized gap will be larger than the gap. Thus the 
Schottky peak moves to the lowest values of the gap in 
respect to the case g s = 0. 



III. NUMERICAL RESULTS 

In recent years, numerical methods have been exten- 
sively developed and applied to quantum many-body 
problems. Most frequently used numerical method for 
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these problems is the exact diagonalization of small sys- 
tems employing the Lanczos technique 2 ^. The exact diag- 
onalization of small correlated systems does not have any 
restrictions on the model. The deficiency of the method 
is in the relative smallness of system sizes. So far the 
method has been essentially restricted to the evaluation 
of the T = static and dynamical quantities, i.e., prop- 
erties of the ground state. 

Jaklic et.al. introduced a method for the evaluation of 
finite-temperature properties, based on the Lanczos diag- 
onalization technique for small systems^. This method, 
is avoid the calculation of all eigenfunctions of the sys- 
tem. Instead, they introduced the procedure where the 
sampling over all states is reduced to a random par- 
tial sampling, while only approximate ground state and 
excited state wave functions, generated by the Lanczos 
technique, are used for the evaluation of matrix elements. 
The size limitations of the method are effectively compa- 
rable to those encountered in the Lanczos-type diago- 
nalization technique applied to the ground state calcula- 
tions. 

In following, we present our numerical results on the 
heat capacity of the several ID spin-1/2 models which 
are obtained by the method of Jaklic. 



A. Alternating Heisenberg Spin-1/2 Chains in a 
Transverse Magnetic Field 

In this section we consider the alternating spin-1/2 
chains in a magnetic field. Since, the antiferromagnctic- 
ferromagnetic (AF-F) chains have a gap in the spin ex- 
citation spectrum, they reveal extremely rich quantum 
behavior in the presence of the magnetic field. 

The ground state phase diagram of the AF-F alter- 
nating chain in a magnetic field is studied by the nu- 
merical diagonalization and the finite-size scaling based 
on the conformal field theory- 2 . It is shown that the 
magnetic state is gapless and described by the Luttinger 
liquid phase. It is also found that the magnetic state 
is characterized by the algebraic decay of the spin cor- 
relation functions. Recently, Yamamoto et.al described 
the magnetic properties of the model in a magnetic field 
in terms of the spinless fermions and the spin waves^. 
They employed the Jordan- Wigner transformation and 
treated the fermionic Hamiltonian within the Hartree- 
Fock approximation. They have also implemented the 
modified spin wave theory to calculate the thermody- 
namic functions as the heat capacity and the magnetic 
susceptibility. More recently, using numerical Lanc- 
zos method, the effect of an uniform transverse magnetic 
field on the ground state phase diagram of a spin-1/2 
AF-F chain with anisotropic ferromagnetic coupling is 
studied 24 . The Hamiltonian of the model under Consid- 
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FIG. 2: a. The excitation gap of a spin-1/2 AF-F chain for 
isotropic case A = 1.0 in the uniform magnetic field for dif- 
ferent size number (TV = 12,16,20). b. The excitation gap 
of a spin-1/2 AF-F chain in the uniform transverse magnetic 
field with anisotropic ferromagnetic coupling A = 0.5 , for 
different size number (N = 12, 16, 20). 

eration on a periodic chain of N sites is given by 

N/2 

H = JaF /"^[ffzj-l^fj + ^2j-1^2j + ^fj-l^fj] 

i=i 

N/2 

- Jf ^ [S%j S^j+l + ^2j^2j+l + ASjy^fj+l] 
N 

+ h^Sf. (7) 

Where Sj' v ' z are spin-1/2 operators on the j-th site. 
Jp and J af denote the ferromagnetic and antiferro- 
magnetic couplings respectively. The limiting case of 
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FIG. 3: a. The Schottky heat capacity peak value, Tm, versus 
the transverse magnetic field, for a spin- 1/2 AF-F chain, and 
isotropic case A = 1.0 in the uniform transverse magnetic 
field for system size: TV = 20. b. The Schottky heat capacity 
peak value, Tm, as a function of the magnetic field, for a 
spin- 1/2 AF-F chain in the uniform transverse magnetic field 
with anisotropic ferromagnetic coupling: A = 0.5, and size 
number: TV = 20. The inset shows the same quantity in the 
intermediate magnetic field region, h ci < h < h C2 . 



isotropic ferromagnetic coupling corresponds to A = 1 
and h is the transverse magnetic field. To explore the 
nature of the excitation spectrum, we use the modified 
Lanczos method to diagonalizc numerically finite chains 
(N = 12,16,20,24). The energies of the few lowest 
eigenstates were obtained for the chains with periodic 
boundary conditions. First, we have computed the three 
lowest energy eigenvalues of TV = 12, 16, 20 chain with 
J af — 2«7p and different values of the anisotropy param- 
eter A. 

In Fig. [2^ we have plotted results of calculations for 
the isotropic case A = 1.0. The excitation gap is 
determined^ in the system as the difference between 
the first excited state and the ground state. As it is 



clearly seen from this figure in the case of zero magnetic 
field the spectrum of the model is gapped. For h ^ 
the gap decreases linearly with h and vanishes at the 
critical field, h Cl — 1.55 ± 0.01. This is the first level 
crossing between the ground state energy and the first 
excited state. To get an accurate estimate of h Cl we 
have obtained the first level crossing for system sizes of 
TV = 12, 14, 24. The finite size behavior of these values 
lead us to h Cl = 1.55 ± 0.01 for TV — ► oo. The spec- 
trum remains gapless for h ci < h < h C2 and becomes 
once again gapped for h > h C2 = 2.0. With increas- 
ing field, for h > h C2 the gap increases linearly with h. 
In the region h Cl < h < h C2 we also observe numerous 
additional level crossing between the lowest eigenstates. 
These level crossing lead to incommensurate effects that 
manifest themselves in the oscillatory behavior of the spin 
correlation functions. All crossings disappear at h > h C2 
and the correlation functions do not contain oscillatory 
terms in this region of the phase diagram. 

In marked contrast with the isotropic case, the similar 
analysis of the few lowest levels for an anisotropic AF-F 
chain in the presence of a transverse magnetic field reveal 
a principally different behavior. The gap as a function 
of the transverse magnetic field h has been computed for 
the anisotropy parameter A = 0.5 and different chain 
lengths TV = 12, 16, 20. In Fig[2b we have plotted results 
of these calculations. As it is seen from the figure, the ex- 
citation spectrum in this case is gapful except at the two 
critical fields h Cl = 1.64 ± 0.01 and h C2 = 1.88 ± 0.01 24 . 
We have employed the phenomenological renormalization 
group (PRG) method^ to determine these critical fields 
(h ci and h C2 ). The PRG equation is 



(TV + 4)gap{N + 4, ti) = Ngap{N, h), 



(8) 



where gap(N, h) = E\ (TV, h) — Eq (TV, h) is the energy gap 
value for chain length TV in a magnetic field h. At the 
critical point, N(E± — Eq) should be size independent 
for large enough systems in which the contribution from 
irrelevant operators is negligible. Thus, we accurately 
determined the critical points by the PRG method. We 
defined h c (N,N + 4) as the TV-dependent fixed point of 
Eq. [51 and it is extrapolated to the thermodynamic limit 
in order to estimate h c . At the critical point h — h! = h c , 
therefore, the curves of N(Ei — Eq) vs h for sizes TV 
and TV + 4 cross at certain values h Cl (TV, TV + 4) and 
h C2 (TV, TV + 4) (finite-size critical points). The thermo- 
dynamic critical points (h ci and h C2 ) are obtained by ap- 
propriately extrapolating h Cl (TV, TV + 4) or h C2 (TV, TV + 4) 
to TV — > oo. 

In the region h Cl < h < h C2 the spin gap, which ap- 
pears at h > h Cl , first increases vs external field and 
after passing a maximum decreases to vanish at h C2 . 
At h > h C2 the gap once again opens and, for a suffi- 
ciently large transverse field becomes proportional to h. 
To study the finite-temperature properties of the model, 
we have used the Jaklic formalism. We have computed a 
hundred lowest eigenvalues of the energies. We have con- 
sidered different values of the transverse magnetic field, 
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FIG. 4: The Schottky heat capacity peak value, Tm, versus 
the energy gap for a spin-1/2 AF-F chain with anisotropic 
ferromagnetic couplings: A = 1.0 and A = 0.5, and size 
number: N = 20. 



h, and anisotropy parameter A. Therefore using these 
hundred eigen-energies, we have computed the heat ca- 
pacity as a function of the temperature (T). The position 
of the Schottky heat capacity peak, Tm is determined as 
Tm- In Fig. [3^ we have plotted Tm as a function of the 
magnetic field h. To arrive at this plot, we have consid- 
ered Jaf = %Jf, A = 1 and N — 20. As it is clearly seen 
from this figure, the position of the Schottky heat capac- 
ity peak, Tm, Tm, decreases by increasing the magnetic 
field up to the critical field, h Cl . In the intermediate 
region of the magnetic filed, h Cl < h < h C2 , Tm is in- 
dependent of magnetic field. The anomaly behavior is 
the result of the finite size effects. With increasing the 
magnetic field, for h > h C2 , Tm increases almost linearly. 

It is surprising which the behavior of Tm versus the 
magnetic field is in complete agreement with the gap be- 
havior respect to the magnetic filed. Thus we conclude 
that the energy gap sign on the position of the Schottky 
heat capacity peak. To confirm our idea, we have plotted 
Tm vs the transverse magnetic field for Jaf = 2Jf, the 
anisotropy parameter A = 0.5 and N = 20 in Fig. [3Jd. 
As we can see from this figure for low transverse mag- 
netic field, h < h Cl , the behavior of Tm is the same 
as the isotropic case. But in the intermediate region, 
h ci < h < h C2 , Tm first increases vs transverse field and 
after passing a maximum, decreases up to h C2 (see the 
inset of the Fig.[3jD). For the higher transverse magnetic 
field, h > h C2 , the position of the Schottky heat capac- 
ity peak increases as a previous one. This behavior of 
the Tm is in complete agreement with the effect of the 
transverse magnetic field on the energy gap (Fig. |5Jd). 

Finally, we have plotted Tm versus the energy gap 
for N — 20, Jaf — 2Jf and different values of the 
anisotropy parameter A = 0.5, 1.0. For convenience, the 



numerical results of the region h < h Cl are showed. It is 
clearly seen, that the possition of the Schottky heat ca- 
pacity peak Tm, increases by increasing the energy gap 
of the system. This is in well agreement with results 
obtained within the two-level model. 



B. The ID AF-Heisenberg model in a staggered 
field 

The general feature developed for the alternating spin- 
1/2 chains in a magnetic field can be applied to the ID 
Heisenberg Hamiltonian with a staggered magnetic field 
h s , 

N 

H = jY^[Si^i+i + h s (-iyS!}. (9) 

i=l 

It is expecte d 26 ' 27 ' 28 ' 29 that the staggered field induces 
an excitation gap in the S = | AF-Heisenberg chain, 
which should be otherwise gapless. The excitation gap 
caused by the staggered field is indeed found in the real 
magnet o 30 ' 31 ! 32 . In the absence of the staggered field 
(h s = 0), the eigenspectra is exactly solvable. In the 
case of the staggered magnetic field (h s ^ 0), the inte- 
grability is lost. The staggered magnetic field produces 
an antiferromagnetic ordered (Neel order) ground-state. 

To examine the effect of the staggered magnetic field on 
the energy gap, we have implemented the modified Lanc- 
zos algorithm for finite-size chains N = 12, 16,20 using 
periodic boundary conditions. The energy gap is deter- 
mined as the difference between the first excited state 
and the ground stated, and calculated for different chain 
lengths and staggered fields h s . The energy gap is deter- 
mined as the difference between the first excited state 
and the ground stated. 

We have plotted, in Fig. [5h., the energy gap versus 
the staggered magnetic field h s . The results have been 
plotted for different chain sizes N = 12, 16, 20. It can be 
seen, that the spectrum is gapless in the absence of the 
staggered magnetic field (h s = 0). The application of a 
staggered magnetic field, induces a gap in the spectrum 
of the model. With increasing the field, for h s > 0, the 
energy gap increases with h s . 

Applying the Jaklic method we have computed a hun- 
dred lowest energies for different values of the staggered 
magnetic field. The heat capacity is computed as a func- 
tion of the temperature. In Fig. [SJa, the position of the 
Schottky anomaly heat capacity peak, Tm is plotted ver- 
sus h s for the chain size N = 20. As it is seen from the 
figure, Tm in this case increases by increasing the stag- 
gered magnetic field. This result is in good agreement 
with the idea that the gap sign on the position of the 
Schottky heat capacity peak, Tm- 
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and numerical approaches. Theoretically, the system is 
mapped to the well known two-level-system (TLS) model. 
In this case, the heat capacity is found exactly as a func- 
tion of the energy gap and the temperature. The position 
of the Schottky heat capacity peak, Tm, is determined. 
It is shown that Tm as a function of the control param- 
eter behaves in the same way as the energy gap versus 
the control parameter. This shows that the gap has an 
influence on the position of the Schottky heat capacity 
peak. 

Numerically, the finite temperature Lanczos method is 
applied. The Lanczos method is implemented to obtain 
a hundred of lowest excited state energies. This formal- 
ism is applied to two model chains up to iV = 24 in 
length. First, the alternating spin- 1/2 chains in a mag- 
netic field are considered. Since, the antiferromagnetic- 
ferromagnetic (AF-F) chains have a gap in the spin exci- 
tation spectrum, they reveal extremely rich quantum be- 
havior in the presence of the magnetic field. The energy 
gap and heat capacity are computed for both isotropic 
and anisotropic cases. The numerical results are com- 
puted for different values of the external magnetic field. 
It is shown, in complete agreement with the theoretical 
results, the field-dependence of the Tm and the energy 
gap are the same. Finally, the ID Heisenberg model with 
a staggered magnetic field h s is investigated. It is shown 
that the staggered field induces an excitation gap in the 
S = 2 AF-Hcisenberg chain, which should be otherwise 
gapless. Using the above numerical procedure, the posi- 
tion of the Schottky heat capacity peak is computed. It 
is confirm, that the energy gap sign on the Schottky heat 
capacity peak. 



FIG. 5: a. The excitation gap of the ID AF-Heisenberg model 
in a staggered field as a function of the magnetic field for 
different size number (N = 12, 16, 20). b. The position of the 
Schottky heat capacity peak, Tm, for the ID AF-Heisenberg 
model in a staggered field versus the staggered magnetic field 
with size number: N = 20. 
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The heat capacity of low-dimensional spin systems is studied using theoretical and numerical 
techniques. Keeping only two energy states, the system is mapped onto the two -level-system (TLS) 
model. Using the low temperature Lanczos method, it is confirmed that the behavior of Tm and the 
energy gap as functions of the control parameter is the same in the two models studied; a conclusion 
that can probably be extrapolated to the general case of any system that possesses an energy gap. 

PACS numbers: 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 

Gaps in the energy spectrum play a crucial role in con- 
densed matter physics. Fundamental properties in su- 
perconductivity or in the fractional quantum Hall effect 
originate from the existence of a gap between the ground 
state and the excited states. In particular, low dimen- 
sional quantum spin systems are extremely interesting to 
study the behavior of the gap. Many exact and numeri- 
cal results on the one dimensional quantum spin systems 
with nearest neighbor couplings have been accumulated 
during last decades. 

The ID spin- 1/2 system has been solved by BetheA 
in 1931 with his famous ansatz. The ansatz allows the 
computation of the energy eigenvalues. The isotropic ID 
spin- 1/2 system with nearest neighbor couplings is gap- 
less. The anisotropic spin- 1/2 chain is denoted by the 
XXZ model 2 -! 3 -. The Hamiltonian of the XXZ model on a 
periodic chain of N sites is 



N 

H = J 2^ SfSf+i 

i=l 



(1) 



where J > is the exchange coupling in the xy easy 
plane, A is the anisotropy in the z direction. The Ising 
regime is governed by A > 1 and there is a gap in the 
excitation spectrum, while for A < — 1, the ground state 
is in the ferromagnetic phase and there is a gap over the 
ferromagnetic state. In the region — 1 < A < 1, the 
ground state of the system, is in the gapless spin-fluid 
phase. 

The spin-1 system is not solvable with the Bcthc ansatz 
or similar techniques. The anisotropic spin-1 chain is 
only gapless 4 - in the region — 1 < A < 0. Haldane^ for- 
mulated in 1983 his famous conjecture that quantum spin 
chains (isotropic) with integer spin S = 1,2,.. have a gap, 
where as chains with half-integer spin S = 1/2, 3/2, ... are 
gapless. 

During the last two decades ladder-systems^ as quan- 
tum spin systems between one and two dimensions have 
been studied as well with numerical methods. Con- 
cerning the ground state properties has been found 
that, ladder-systems with an even number of legs (I = 



2,4,6,...) have a gap; those with an odd number 
(I = 1,3,5,...) do not. In particular, since the an- 
tiferromagnetic two-leg ladder systems have a gap in 
the spin excitation spectrum, they reveal extremely 
rich quantum behavior in the presence of a magnetic 
field^. Such quantum phase transitions in spin systems 
with gapped excitation spectrum were indeed studied 
experimentall y 8 ! 9 ! 10 ! 11 ! 12 ! 13 . 

On the other hand, investigating the behavior of the 
energy gap of spin systems in vicinity of quantum critical 
points has attracted much interest recentl y 3 ' 14 ! 15 ' 16 . In 
general, the critical point of an thermodynamic system 
in the Hamiltonian formulation is defined as the value at 
which the energy gap vanishes as a power low, which is 
known as the scaling behavior. The opening of the energy 
gap in the vicinity of the quantum critical point is found 
to scale with a critical exponent. The value of the critical 
gap exponent is very important to find the universality 
class of a continuous quantum phase transition. 

The discovery of gapless or gaped excitations have 
led to the investigation of the thermodynamic proper- 
ties. One of the most important thermodynamic func- 
tions is the heat capacity. Usually, there is a lambda-type 
anomaly in figure of the heat capacity versus tempera- 
ture. It was interpreted as indicating a phase transition 
to a magnetically ordered phase. An important charac- 
teristic of the low-dimensional magnets is the absence of 
the long range order in models with a continuous sym- 
metry at any finite temperature 1 ^. There is also a broad 
maximum in the plot of the heat capacity vs tempera- 
ture, characteristic of low dimensional systems. This is 
known as the Schottky peak. 

In this paper, theoretical and numerical results are re- 
ported for the low-temperature behavior of the heat ca- 
pacity in low dimensional spin systems. Theoretically, 
by keeping only two lowest states, the system is mapped 
to the well known two-level-system (TLS) model. In this 
case, the heat capacity is found exactly as a function 
of the energy gap and the temperature. It is shown that 
the position of the Schottky heat capacity peak, Tm , and 
the energy gap behaves in the same way as a function of 
the control parameter. In Section II, the mapping to the 
TLS model is explained and the theoretical results are 
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presented. In Section III, the results of the low temper- 
ature Lanczos calculations are presented. Numerically, 
Tm as a function of the control parameter is computed 
for "the alternating spin- 1/2 chains in a magnetic field 
h" and "the ID Heisenberg Hamiltonian with a staggered 
magnetic field h s " ■ Finally, the summary and conclusions 
are presented in Section IV 



II. LOW TEMPERATURE LIMIT HEAT 
CAPACITY: TWO LEVEL SYSTEM APPROACH 

In this section we discuss a theoretical approach to find 
the effect of the energy gap on the heat capacity of the 
quantum spin systems. Spectrum energy of a quantum 
system may be gapful or gapless. Since, we are going to 
study the sign of the gap on the heat capacity, gapful 
systems are considered. At very low temperature we can 
consider only lowest energy levels. If we keep only two 
lowest energy states, the system maps to the two level 
system (TLS) modeli&. We have assumed that the energy 
gap of this TLS is g where E\ = Eq + g, E and E\ are 
the ground and first exited states respectively. 

Our purpose is to determine the behavior of the heat 
capacity for the quantum spin systems. The heat capac- 
ity is expressed by the following relation 



C v — 



1 



K B T 2 



d 2 InZ 



dp 2 



where Z is the partition function denoted by 
Z = Trje-^}, 



(2) 



(3) 



Here (3 = K T , Kb is the Boltzmann constant and T is 
temperature. We have assumed that Kb = 1. At very 
low temperature limit (TLS limit), we can write 
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And finally from the above equation one can shows 

,.2 



Cy 



cosh 2 (x) ' 



(4) 



(5) 



where we defined x 



2T • 



This result predict a Schot- 



tky like peak of the heat capacity behavior versus the x 
variation. The position of the Schottky peak takes place 
at xm — 1-2, where we have: tanh(a;jvf) = 1/- t m- This 
result shows explicitly an upward increase of the heat 
capacity versus the magnetic field for x < xm and mono- 
tonic decrease for x > xm- We have plotted the thermal 
behavior of Eq. [5] in Fig. [TJ Here, xm corresponds to 
the position of the Schottky peak, Tm, in a constant gap 
value. It is clear that by increasing the gap value, \g\, 
the peak approach to higher temperature and inversely. 
Therefore we can conclude that width of the energy gap 
may affect the position of the Schottky heat capacity 
peak. At very low temperature regime, all gapped ther- 
modynamic systems (N — > oo, N = number of the spins) 
can be mapped to the above TLS model. 




Temperature 

FIG. 1: Temperature dependence of the heat capacity of TLS 
(Eq. [SJ) for different energy gap values (g). Tm shows the 
position of the Schottky heat capacity peak for the system 
with g = 2. 



Up to now we did not consider any degeneracy. In the 
general case both first two energy levels (TLS) have a 
degeneracy. It is well known for two level system (TLS) 
that the account of degeneracy leads to change the effec- 
tive gap and consequently to change the Tm for Schottky 
anomaly of heat capacity 19 . In peresent of degeneracy by 
a little of manipulation, one can shows 



Cy — 



cosh 2 (a; + x') ' 



(6) 



which is very similar to Eq. [5J Here we denoted x' — — 



IT ' 



and g s = T In ^ . Where d\ and di have been con- 
sidered for the order of the degeneracy of the ground 
state and first exited state of the system respectively. 
As we mentioned it before this result predict a Schot- 
tky like peak in the heat capacity behavior. The Schot- 
tky heat capacity peak takes place at x m which satisfied 
tanh(a;M + x' M ) = 1/xm, where x' M — 27^7- If we as ~ 
sumed that d± < di then g a < and therefore the gener- 
alized gap (g + g s ) will be smaller than the gap. Thus the 
Schottky peak moves to the higher values of the gap in 
respect to the case g s — (without degeneracy) . In same 
manner for the case: di < d\ we have g s > 0, therefore 
the generalized gap will be larger than the gap. Thus the 
Schottky peak moves to the lowest values of the gap in 
respect to the case g s = 0. 



III. NUMERICAL RESULTS 

In recent years, numerical methods have been exten- 
sively developed and applied to quantum many-body 
problems. Most frequently used numerical method for 
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these problems is the exact diagonalization of small sys- 
tems employing the Lanczos technique 2 ^. The exact diag- 
onalization of small correlated systems does not have any 
restrictions on the model. The deficiency of the method 
is in the relative smallness of system sizes. So far the 
method has been essentially restricted to the evaluation 
of the T = static and dynamical quantities, i.e., prop- 
erties of the ground state. 

Jaklic et.al. introduced a method for the evaluation of 
finite-temperature properties, based on the Lanczos diag- 
onalization technique for small systems^. This method, 
is avoid the calculation of all eigenfunctions of the sys- 
tem. Instead, they introduced the procedure where the 
sampling over all states is reduced to a random par- 
tial sampling, while only approximate ground state and 
excited state wave functions, generated by the Lanczos 
technique, are used for the evaluation of matrix elements. 
The size limitations of the method are effectively compa- 
rable to those encountered in the Lanczos-type diago- 
nalization technique applied to the ground state calcula- 
tions. 

In following, we present our numerical results on the 
heat capacity of the several ID spin-1/2 models which 
are obtained by the method of Jaklic. 



A. Alternating Heisenberg Spin-1/2 Chains in a 
Transverse Magnetic Field 

In this section we consider the alternating spin-1/2 
chains in a magnetic field. Since, the antiferromagnctic- 
ferromagnetic (AF-F) chains have a gap in the spin ex- 
citation spectrum, they reveal extremely rich quantum 
behavior in the presence of the magnetic field. 

The ground state phase diagram of the AF-F alter- 
nating chain in a magnetic field is studied by the nu- 
merical diagonalization and the finite-size scaling based 
on the conformal field theory- 2 . It is shown that the 
magnetic state is gapless and described by the Luttinger 
liquid phase. It is also found that the magnetic state 
is characterized by the algebraic decay of the spin cor- 
relation functions. Recently, Yamamoto et.al described 
the magnetic properties of the model in a magnetic field 
in terms of the spinless fermions and the spin waves^. 
They employed the Jordan- Wigner transformation and 
treated the fermionic Hamiltonian within the Hartree- 
Fock approximation. They have also implemented the 
modified spin wave theory to calculate the thermody- 
namic functions as the heat capacity and the magnetic 
susceptibility. More recently, using numerical Lanc- 
zos method, the effect of an uniform transverse magnetic 
field on the ground state phase diagram of a spin-1/2 
AF-F chain with anisotropic ferromagnetic coupling is 
studied 24 . The Hamiltonian of the model under Consid- 
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FIG. 2: a. The excitation gap of a spin-1/2 AF-F chain for 
isotropic case A = 1.0 in the uniform magnetic field for dif- 
ferent size number (TV = 12,16,20). b. The excitation gap 
of a spin-1/2 AF-F chain in the uniform transverse magnetic 
field with anisotropic ferromagnetic coupling A = 0.5 , for 
different size number (N = 12, 16, 20). 

eration on a periodic chain of N sites is given by 

N/2 

H = JaF /"^[ffzj-l^fj + ^2j-1^2j + ^fj-l^fj] 

i=i 

N/2 

- Jf ^ [S%j S^j+l + ^2j^2j+l + ASjy^fj+l] 
N 

+ h^Sf. (7) 

Where Sj' v ' z are spin-1/2 operators on the j-th site. 
Jp and J af denote the ferromagnetic and antiferro- 
magnetic couplings respectively. The limiting case of 



4 



(al 



N = 2 
J AK = 2J r 
A=l .0 



_i i i i I i i i i I i i i i L_ 



(b) 




N = 20 

A = 0.5 

O 

O 



^ooo 



_1 I I I I I I I I I I I I I l_ 



FIG. 3: a. The Schottky heat capacity peak value, Tm, versus 
the transverse magnetic field, for a spin- 1/2 AF-F chain, and 
isotropic case A = 1.0 in the uniform transverse magnetic 
field for system size: TV = 20. b. The Schottky heat capacity 
peak value, Tm, as a function of the magnetic field, for a 
spin- 1/2 AF-F chain in the uniform transverse magnetic field 
with anisotropic ferromagnetic coupling: A = 0.5, and size 
number: TV = 20. The inset shows the same quantity in the 
intermediate magnetic field region, h ci < h < h C2 . 



isotropic ferromagnetic coupling corresponds to A = 1 
and h is the transverse magnetic field. To explore the 
nature of the excitation spectrum, we use the modified 
Lanczos method to diagonalizc numerically finite chains 
(N = 12,16,20,24). The energies of the few lowest 
eigenstates were obtained for the chains with periodic 
boundary conditions. First, we have computed the three 
lowest energy eigenvalues of TV = 12, 16, 20 chain with 
J af — 2«7p and different values of the anisotropy param- 
eter A. 

In Fig. [2^ we have plotted results of calculations for 
the isotropic case A = 1.0. The excitation gap is 
determined^ in the system as the difference between 
the first excited state and the ground state. As it is 



clearly seen from this figure in the case of zero magnetic 
field the spectrum of the model is gapped. For h ^ 
the gap decreases linearly with h and vanishes at the 
critical field, h Cl — 1.55 ± 0.01. This is the first level 
crossing between the ground state energy and the first 
excited state. To get an accurate estimate of h Cl we 
have obtained the first level crossing for system sizes of 
TV = 12, 14, 24. The finite size behavior of these values 
lead us to h Cl = 1.55 ± 0.01 for TV — ► oo. The spec- 
trum remains gapless for h ci < h < h C2 and becomes 
once again gapped for h > h C2 = 2.0. With increas- 
ing field, for h > h C2 the gap increases linearly with h. 
In the region h Cl < h < h C2 we also observe numerous 
additional level crossing between the lowest eigenstates. 
These level crossing lead to incommensurate effects that 
manifest themselves in the oscillatory behavior of the spin 
correlation functions. All crossings disappear at h > h C2 
and the correlation functions do not contain oscillatory 
terms in this region of the phase diagram. 

In marked contrast with the isotropic case, the similar 
analysis of the few lowest levels for an anisotropic AF-F 
chain in the presence of a transverse magnetic field reveal 
a principally different behavior. The gap as a function 
of the transverse magnetic field h has been computed for 
the anisotropy parameter A = 0.5 and different chain 
lengths TV = 12, 16, 20. In Fig[2b we have plotted results 
of these calculations. As it is seen from the figure, the ex- 
citation spectrum in this case is gapful except at the two 
critical fields h Cl = 1.64 ± 0.01 and h C2 = 1.88 ± 0.01 24 . 
We have employed the phenomenological renormalization 
group (PRG) method^ to determine these critical fields 
(h ci and h C2 ). The PRG equation is 



(TV + 4)gap{N + 4, ti) = Ngap{N, h), 



(8) 



where gap(N, h) = E\ (TV, h) — Eq (TV, h) is the energy gap 
value for chain length TV in a magnetic field h. At the 
critical point, N(E± — Eq) should be size independent 
for large enough systems in which the contribution from 
irrelevant operators is negligible. Thus, we accurately 
determined the critical points by the PRG method. We 
defined h c (N,N + 4) as the TV-dependent fixed point of 
Eq. [51 and it is extrapolated to the thermodynamic limit 
in order to estimate h c . At the critical point h — h! = h c , 
therefore, the curves of N(Ei — Eq) vs h for sizes TV 
and TV + 4 cross at certain values h Cl (TV, TV + 4) and 
h C2 (TV, TV + 4) (finite-size critical points). The thermo- 
dynamic critical points (h ci and h C2 ) are obtained by ap- 
propriately extrapolating h Cl (TV, TV + 4) or h C2 (TV, TV + 4) 
to TV — > oo. 

In the region h Cl < h < h C2 the spin gap, which ap- 
pears at h > h Cl , first increases vs external field and 
after passing a maximum decreases to vanish at h C2 . 
At h > h C2 the gap once again opens and, for a suffi- 
ciently large transverse field becomes proportional to h. 
To study the finite-temperature properties of the model, 
we have used the Jaklic formalism. We have computed a 
hundred lowest eigenvalues of the energies. We have con- 
sidered different values of the transverse magnetic field, 
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FIG. 4: The Schottky heat capacity peak value, Tm, versus 
the energy gap for a spin-1/2 AF-F chain with anisotropic 
ferromagnetic couplings: A = 1.0 and A = 0.5, and size 
number: N = 20. 



h, and anisotropy parameter A. Therefore using these 
hundred eigen-energies, we have computed the heat ca- 
pacity as a function of the temperature (T). The position 
of the Schottky heat capacity peak, Tm is determined as 
Tm- In Fig. [3^ we have plotted Tm as a function of the 
magnetic field h. To arrive at this plot, we have consid- 
ered Jaf = %Jf, A = 1 and N — 20. As it is clearly seen 
from this figure, the position of the Schottky heat capac- 
ity peak, Tm, Tm, decreases by increasing the magnetic 
field up to the critical field, h Cl . In the intermediate 
region of the magnetic filed, h Cl < h < h C2 , Tm is in- 
dependent of magnetic field. The anomaly behavior is 
the result of the finite size effects. With increasing the 
magnetic field, for h > h C2 , Tm increases almost linearly. 

It is surprising which the behavior of Tm versus the 
magnetic field is in complete agreement with the gap be- 
havior respect to the magnetic filed. Thus we conclude 
that the energy gap sign on the position of the Schottky 
heat capacity peak. To confirm our idea, we have plotted 
Tm vs the transverse magnetic field for Jaf = 2Jf, the 
anisotropy parameter A = 0.5 and N = 20 in Fig. [3Jd. 
As we can see from this figure for low transverse mag- 
netic field, h < h Cl , the behavior of Tm is the same 
as the isotropic case. But in the intermediate region, 
h ci < h < h C2 , Tm first increases vs transverse field and 
after passing a maximum, decreases up to h C2 (see the 
inset of the Fig.[3jD). For the higher transverse magnetic 
field, h > h C2 , the position of the Schottky heat capac- 
ity peak increases as a previous one. This behavior of 
the Tm is in complete agreement with the effect of the 
transverse magnetic field on the energy gap (Fig. |5Jd). 

Finally, we have plotted Tm versus the energy gap 
for N — 20, Jaf — 2Jf and different values of the 
anisotropy parameter A = 0.5, 1.0. For convenience, the 



numerical results of the region h < h Cl are showed. It is 
clearly seen, that the possition of the Schottky heat ca- 
pacity peak Tm, increases by increasing the energy gap 
of the system. This is in well agreement with results 
obtained within the two-level model. 



B. The ID AF-Heisenberg model in a staggered 
field 

The general feature developed for the alternating spin- 
1/2 chains in a magnetic field can be applied to the ID 
Heisenberg Hamiltonian with a staggered magnetic field 
h s , 

N 

H = jY^[Si^i+i + h s (-iyS!}. (9) 

i=l 

It is expecte d 26 ' 27 ' 28 ' 29 that the staggered field induces 
an excitation gap in the S = | AF-Heisenberg chain, 
which should be otherwise gapless. The excitation gap 
caused by the staggered field is indeed found in the real 
magnet o 30 ' 31 ! 32 . In the absence of the staggered field 
(h s = 0), the eigenspectra is exactly solvable. In the 
case of the staggered magnetic field (h s ^ 0), the inte- 
grability is lost. The staggered magnetic field produces 
an antiferromagnetic ordered (Neel order) ground-state. 

To examine the effect of the staggered magnetic field on 
the energy gap, we have implemented the modified Lanc- 
zos algorithm for finite-size chains N = 12, 16,20 using 
periodic boundary conditions. The energy gap is deter- 
mined as the difference between the first excited state 
and the ground stated, and calculated for different chain 
lengths and staggered fields h s . The energy gap is deter- 
mined as the difference between the first excited state 
and the ground stated. 

We have plotted, in Fig. [5h., the energy gap versus 
the staggered magnetic field h s . The results have been 
plotted for different chain sizes N = 12, 16, 20. It can be 
seen, that the spectrum is gapless in the absence of the 
staggered magnetic field (h s = 0). The application of a 
staggered magnetic field, induces a gap in the spectrum 
of the model. With increasing the field, for h s > 0, the 
energy gap increases with h s . 

Applying the Jaklic method we have computed a hun- 
dred lowest energies for different values of the staggered 
magnetic field. The heat capacity is computed as a func- 
tion of the temperature. In Fig. [SJa, the position of the 
Schottky anomaly heat capacity peak, Tm is plotted ver- 
sus h s for the chain size N = 20. As it is seen from the 
figure, Tm in this case increases by increasing the stag- 
gered magnetic field. This result is in good agreement 
with the idea that the gap sign on the position of the 
Schottky heat capacity peak, Tm- 
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and numerical approaches. Theoretically, the system is 
mapped to the well known two-level-system (TLS) model. 
In this case, the heat capacity is found exactly as a func- 
tion of the energy gap and the temperature. The position 
of the Schottky heat capacity peak, Tm, is determined. 
It is shown that Tm as a function of the control param- 
eter behaves in the same way as the energy gap versus 
the control parameter. This shows that the gap has an 
influence on the position of the Schottky heat capacity 
peak. 

Numerically, the finite temperature Lanczos method is 
applied. The Lanczos method is implemented to obtain 
a hundred of lowest excited state energies. This formal- 
ism is applied to two model chains up to iV = 24 in 
length. First, the alternating spin- 1/2 chains in a mag- 
netic field are considered. Since, the antiferromagnetic- 
ferromagnetic (AF-F) chains have a gap in the spin exci- 
tation spectrum, they reveal extremely rich quantum be- 
havior in the presence of the magnetic field. The energy 
gap and heat capacity are computed for both isotropic 
and anisotropic cases. The numerical results are com- 
puted for different values of the external magnetic field. 
It is shown, in complete agreement with the theoretical 
results, the field-dependence of the Tm and the energy 
gap are the same. Finally, the ID Heisenberg model with 
a staggered magnetic field h s is investigated. It is shown 
that the staggered field induces an excitation gap in the 
S = 2 AF-Hcisenberg chain, which should be otherwise 
gapless. Using the above numerical procedure, the posi- 
tion of the Schottky heat capacity peak is computed. It 
is confirm, that the energy gap sign on the Schottky heat 
capacity peak. 



FIG. 5: a. The excitation gap of the ID AF-Heisenberg model 
in a staggered field as a function of the magnetic field for 
different size number (N = 12, 16, 20). b. The position of the 
Schottky heat capacity peak, Tm, for the ID AF-Heisenberg 
model in a staggered field versus the staggered magnetic field 
with size number: N = 20. 
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